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We systematically investigate the universal spin decoherence dynamics of a localized electron 
in an arbitrary nuclear spin bath, which can be even far away from equilibrium due to the weak 
nuclear-lattice interaction. We show that the electron spin relaxation dynamics (as well as spin pure 
dephasing and Hahn echo decay) can always have a universal behavior as long as the initial state is 
composed of a sufficiently large amount of spin eigenstates. For a given system, the pattern of the 
universal dynamics depends on the complicated initial condition only via a single parameter, which 
measures the amount of phase coherence between different spin eigenstates in the initial state. Our 
results apply even when the number of the involved nuclei is not large, and therefore provide a solid 
foundation in the comparison of theoretical/numerical results with the experimental measurement. 
As an example, we also show numerical results for systems of noninteracting spin bath in zero 
magnetic held regime, and discuss the features of universal decoherent dynamics. 



I. INTRODUCTION 

Localized spins in solid-state systems are one of the 
most promising candidates for realizing the quantum 
computation due to its long coherence time [l[ and the 
possible scalability |2| . Recently, quantum control of sin- 
gle localized spin becomes experimentally feasible Q, but 
the nuclear spin bath induced decoherence, which is the 
dominant decoherent mechanism in the low temperature 
regime, is still hindering the further developments. To 
study the effects of nuclear spin bath, both analytical ap- 
proaches [100] and numerical simulations are devel- 
oped for different parameter regime. The effect of dipolar 
interaction between nuclear spins are also studied in Refs. 
@. However, to the best of our knowledge, an uncon- 
tested conclusion about the spin decoherence dynamics 
and its relation to the experimental measurement is still 
unavailable, even though the deleterious effects of nuclear 
spin have been verified in recent experiments [lOMllj . 

From experimental side, the most crucial limitation 
results from the fact that the initial nuclear spin con- 
figuration is very little known nor controllable. This is 
a highly nontrivial problem, because even if a thermal- 
ized spin bath is assumed in the beginning (as done in 
most theoretical work [5.6,9]), any quantum measure- 
ment or manipulation of electron spin can just destroy 
the equilibrium and lead to a highly non-equilibrium nu- 
clear spin dynamics. The coherent time of nuclear spin 
bath is known to be extremely long (can exceed Is in 
GaAs quantum well [l2j and 25s in 29 Si isotope [Hj]) and 
therefore it is very questionable if the nuclear spin bath 
could be wcll-thermalized for the next quantum measure- 
ment/manipulation in a short time during the quantum 
computation process. In order to have a meaningful com- 
parison between theoretical results and the experimental 
measurement, the first and the most important question 
one should ask is if there could be any universal dynam- 
ics in such a system, which is insensitive to the details of 
initial nuclear spin configuration. 

From theoretical side, answering above question is also 
very difficult because the spin dynamics of one configura- 



tion can be very different from the other Q even though 
their initial configurations are similar. Moreover, in a 
typical quantum dot system, the number of nuclei can be 
very huge (N ~ 10 3-5 ), and hence it is also a significant 
challenge for ordinary numerical simulation to explore 
such huge phase space. These challenges are fundamen- 
tally important to the understanding of the spin decoher- 
ence mechanism and to its future application in quantum 
computation. However, to the best of our knowledge, 
there has no systematically study in the literature to this 
important issue. 

In this paper we address this issue by rigorously prove 
the existence of a generic and universal electron spin de- 
coherent dynamics in an arbitrary nuclear spin bath. By 
"universal dynamics" we mean an electron spin evolution 
which is of zero standard deviation over different initial 
condition in the whole phase space. More precisely, we 
show that (1) the universality of spin decoherence always 
exists if only the initial state is composed of sufficient 
large amount of spin eigenstates, and (2) for a given sys- 
tem, such universal dynamics depends on the initial con- 
figuration only through a single parameter, which mea- 
sures the amount of phase coherence between spin eigen- 
states of the initial wavefunction. (3) The universality 
is ensured by the large amount of phase space rather 
than the large value of nuclear number, N, and there- 
fore numerical simulation for a small size system (say 
N ~ 1 — 20) can be still good enough to compare with a 
realistic system of much more nuclei [l4{ . Finally, (4) the 
universality of spin dynamics applies to the decoherence 
of the diagonal part (S z ) as well as the off-diagonal part 
{S x ) of electron spin, no matter in a free induction decay 
(FID) or in a Hahn echo decay. Therefore our results 
resolves the fundamental problems in the comparison of 
a theoretical calculation and an experimental measure- 
ment, and provide a new direction for the future study 
of the spin decoherence. We also study the spin dynamics 
for systems of different electron/nuclear spins, and find 
that the spin dynamics is mainly determined by the ge- 
ometric structure of the system density of states and is 
therefore insensitive to the magnitude of nuclear spin. 



This paper is organized as fellows: In section II we de- 
scribe the system Hamiltonian and the initial wavefunc- 
tions in our study. In section III we show the universal 
dynamics of electron spin relaxation. We study the uni- 
versal dynamics by using both numerical and analytical 
methods. In section IV we discuss the microscopic origin 
of the universality. In section V, we generalize our con- 
sideration to other spin systems. We conclude in section 
VI. 



II. SPIN EIGENSTATES AND PHASE SPACE 

A general spin decoherence due to nuclear spin bath is 
described by the following Hamiltonian: 

N 

H = S-J2 A ^i + Hn-n + H Z , (1) 

i 

where S and 1^ are respectively the dimensionless spin op- 
erators (h = 1) of the localized electron and the nucleus 
at lattice site i. At is the hyperfine coupling strength, de- 
pending on the wavefunction profile of the localized elec- 
tron (Fig. HJa)), and we use Ai = Aoe - ^ 1 ^ for our 
numerical calculation with N being the number of nu- 
clei. We note that changing N will not change the 
shape of the electron wavefunction, and therefore 
increasing number of nuclear just reduce the stan- 
dard deviation of the spin dynamics instead of its 
average value. H n -n and Hz are respectively the inter- 
action between nuclear spin and the Zeeman term due to 
external magnetic field. Even in the simplest case, where 
both Hz and H n - n are zero or neglected, the result- 
ing dynamics due to the electron-nuclear coupling only is 
still quite complex, because it involves a huge amount of 
eigenstates in the Hilbert space. In order to have a mean- 
ingful comparison between theoretical(numerical) results 
and the experimental observation, the first question one 
should ask is if there could be any universal dynamics in 
such spin system, which is insensitive to a general initial 
condition of the system and therefore can be observed 
and repeatable in a realistic experiment. After all, it is 
very difficult to control and/or manipulate the spin con- 
figuration of nuclei in solid state systems. This question, 
to the best of our knowledge, has not been answered or 
even not addressed yet in the literature. 

It is convenient to use spin eigenstate, |S) e ® \j) n = 
|S) e CB5 |{ii, z ,i2,zj " ■ ,lN,z}j)n, as the basis of calcula- 
tion, where |S) e is electron spin eigenstate along certain 
direction (will be specified below) and Ji iZ is the nuclear 
spin eigenvalue along the direction of magnetic field (£) 
at the ith site. For simplicity, in this paper, we assume 
the electron spin is initially polarized only along z or x 
axis, and therefore a general initial wavefunction can be 
written to be: \ipo) x , z = \S x>z ) e <g> Y,j=i a j\j)n, where 
a j = r j£ l(pj is the coefficient of the jth spin eigenstate 
with phase ifj and amplitude rj. Here \S z ) e = |+} e 
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FIG. 1: (a) Schematic pictures for electron spin (white cir- 
cle) coupled to a nuclear spin bath (black dots). The elec- 
tron is assumed to be described by the orbital envelope wave 
function cf>, and interacts with the nuclear spins (located at 
Pi) via a hyperfine interaction At — Ao ^(r^) | 2 where Aq is 
the coupling strength, (b) and (c) are two spin eigenstates 
with maximum/minimum energies in zero magnetic field and 
J z = case. 

and \S x ) e = (|+) e + |— }e), and is the size of the 
Hilbert space of nuclear bath. Using above expression, 
we consider the following three spin dynamics, which 
are related to the spin relaxation, spin pure dephasing, 
and Hahn echo decay respectively. The first two can 
be expressed as (S ZjX (t)) = z ,x(^o\S z , x (t)\ip ) ZtX = 
Y^r^S^it) + E^r 3 r ie -^~^S;f(t), where 
S*f(t) = n (j\ <S) e{S z , x \S z , x {t)\S z . x ) e <g> \l) n is the matrix 
element. (S z (t)) can be very different from (S x (t)) if 
the nuclear spin is polarized by external magnetic field 
or with finite total angular momentum in a certain 
direction. Similarly the Hahn echo decay is given by 
p+_(r) = (,{+\ph{t)\— ) e , where the Hahn echo den- 
sity matrix p H (r) = Tr n {U(T)\ipo) xx (ip \U{T^} = 
T, jn {j\U{r)\^),x{MU{ry\3)n and U(t) 
e -iHr G ^ e -iHT |g|[x5| characteristic time scale 

T 2 of pure dephasing is related to the single spin FID, 
while Hahn echo decay 15] is usually used to extract 
single spin behavior from an ensemble measurement. 

III. UNIVERSAL DYNAMICS 

In this section we show that the universality of spin 
decoherence always exists if only the initial state is com- 
posed of sufficient large amount of spin eigenstates, and 
for a given system, such universal dynamics depends on 
the initial configuration only through a single parame- 
ter, which measures the amount of phase coherence be- 
tween spin eigenstates of the initial wavefunction. Wc 
first show the numerical results for spin relaxation, spin 
pure dephasing, and Hahn echo decay respectively. Then 
we rigorously give the proof of the universality. 

A. Numerical study 

In order to explore the spin dynamics from different 
initial conditions in the whole phase space, in this pa- 
per we allow both the amplitude, {V-,-}, and the phase, 
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{<Pj}, to be independent variables and randomly chosen 
according to distribution functions P r { r j) an d V v ((pj) 
respectively. The ensemble-averaged spin dynamics for 
(S z (t)) becomes [(S z (t))] = [fe$fe > where [f(r)] r = 
f Q V r {r)f{r)dr denotes the average of a function f(r), 

and similarly [/(</?)] v ee J Q 27r V v {(p)f{ip)dip. [{S z (Q))] r>lp in 
the denominator is for normalization. At the same time, 
the associated normalized standard deviation (NSD) is 
defined as follows: 



a(t) 



[{S z (t)V 2 



(2) 



1',-p 



Similar definition of averaged dynamics as well as the 
NSD for (S x (t)) and p+_(t) can be obtained easily. We 
note that, if the NSD of the averaged spin dynamics goes 
to zero in the limit of infinite phase space, the averaged 
dynamics is also "the most probable" dynamics with al- 
most zero probability in the other time-evolution behav- 
ior. As a result, we can define it as a universal dynam- 
ics of the given system, independent (in the probability 
sense) of the details of initial nuclear spin configuration. 
On the other hand, the system has no universal dynam- 
ics if the NSD is of the order of one, since the average 
value could not represent the characteristic dynamics of 
a general initial condition. 

Before analytically studying the universality of spin 
dynamics in a general system, it is more instructive 
to show some numerical results of the simplest sys- 
tem without magnetic field and nuclear spin interaction 
(Hz = H n - n = 0). We will first present the result 
for (S z (t)) then the results for (S x (t)) and pf_(t). For 
the convenience of later discussion, we restrict the cal- 
culation inside a subspace, T, where the total angular 
momentum, J z — S z + is zero, and choose 

V r (r) = 7 +(l- 7 )5(r) and 7^) = (1-£)/2tt + £%). 
Here 7 S [0, 1] can be understood as the probability to 
have a nonzero contribution in the subspace, T, while 
£ £ [0, 1] is the probability to have a phase coherence at 
a given value (set to be zero). They satisfys the normal- 
ization condition: J V r {r)dr = J Q V P v (if)d(p = 1 for all 
£ and 7. 

In Fig. Hta)-(c) we show the averaged electron spin 
(S = |) relaxation [(S* (*))], in a noninteracting spin 
bath (7 = ~) with zero magnetic field for £ = 0, 0.5, and 
1 with two different values of 7. We observe that when 
7 is small (7 = 0.016), meaning only a few spin eigen- 
states are involved in \ipo) z , the NSDs are very large, 
i.e., no universal dynamics. This explains why in the 
literature different initial states can result in very differ- 
ent time-evolution patterns pj. When 7 becomes larger 
(7 = 0.32), the NSD decreases in all figures (see also 
Fig. [2Jd)), while the averaged dynamics also begin to 
show different behavior for different value of £: A two- 
decay curves for £ — » and a single mode oscillation for 
£ — > 1. In fact, as we will show later, the NSD always 
decrease to zero even 7 is finite, as long as the size of 
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FIG. 2: (Color online) (a)-(c) Averaged electron spin re- 
laxation (symboled lines) for £=0, 0.5, and 1. respectively. 
Dashed and dotted lines are the uncertain range ([(S'(t))] ± 
a(t)) for 7 = 0.016 and 0.32 respectively, (d) Time averaged 
NSD (a = limT-,00 a(t)dt/T) v.s. 7. N = 11 in all figures. 



phase space becomes large enough. These results indi- 
cate that a universal dynamics can always be expected if 
the initial state is composed of a sufficiently large portion 
of spin eigenstates in the phase space, simply due to the 
strong quantum interference effects. We could also show 
that the two decay time scale of Fig. 2(a) is due the the 
structure of the system density of states, and wil discuss 
that in more details in the latter section. 

In Fig. [3fa) and (b) we respectively plot pure dephas- 
ing ((S x (t)}) and Hahn echo decay (/9+_(t)) as a function 
of time for fully coherent (£ = 1, filled circle) and fully 
incoherent (£ = 0, open circle) choice of initial condi- 
tion. To simplify the numerical calculation, we choose 
an initial wavefunction in a subspace of X^ili h,* = 
for the pure dephasing ((S x (t)}, (a) ), and J z — (same 
as (S z (t))) for the Hahn echo decay (p+_(r), (b)). The 
electron for the former case is initially polarized in the x 
direction so that the average total angular momentum in 
z direction, (S z (t)), is still zero. As a result, the dynam- 
ics of dephasing, (S x (t)), is different from the relaxation, 
(S z (t)) due to the different choice of subspace where the 
initial wavefunction is defined. We believe such convec- 
tion is justified and will not affect any of our conclusion, 
because here we just use this numerical results as an ex- 
ample to understand the general properties of the uni- 
versal dynamics. Full numerical results for any realistic 
situation will need a much larger phase space and much 
longer time. Within this subspace, different initial wave- 
functions still result in different dynamics ( not show here 
). However, when the initial wavefunction is composed of 
sufficient large amount of eigenstates in the subspace, we 
again find a universal dynamics with almost zero NSD. 
In Fig 3(a) and (b), we show results for 7 = 1 for pure 
dephasing and Hahn echo decay in the two subspace de- 
scribed above. In (c), the NSD of the Hahn echo decay 



4 



(Sd(t)) shows a clear two-decay process: with a fast de- 
cay in short time and a slow decay in long time. How- 
ever, the off-diagonal part (S (t)) does not decay at all, 
and shows a single mode oscillation. It is easy to see 
that the numerical results shown in Fig. [2ja)-(c) can be 
obtained as a superposition of Sd{t) and S (t), just as 
suggested by Eq. ([3]). Numerical comparison between 
these two approaches (ensemble average before and after 
time-revolution) agree excellently well (not shown here), 
showing that only a single parameter, (5, is necessary to 
reproduce all the ensemble averaged spin relaxation dy- 
namics. 

In order to exam if the ensemble-averaged results of 
Eq. © is a universal dynamics, we need to calcu- 
late the fluctuation (NSD, Eq. ©) about this aver- 
age. For simplicity, we first study the case in the com- 
FIG. 3: (Color online) (a) [(S x (t))] and (b) p^_(r) as a func- p l e tely random phase limit, say [e l *] v = 0. We then have 
tion of time t and r respectively (see text). Lines with filled 
symbols and open symbols are for £ = 1 and £ = 0. (c) The 
NSD at t — 10 for Hahn echo decay, (d) Sj(t) (open circle) 
and S'o(t) (filled circle) for the spin relaxation dynamics in 
Eq. ©. 
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is plotted as function of 7 at r = 10. From these results, 
we find similar single mode oscillation for £ = 1, while 
a two-decay curves for £ = in all the three dynamics 
((S z (t)),(S x (t)) and pf_(T)). 



B. Analytical study 

To analytically study the universal spin dynamics, we 
have to do the ensemble-average first so that 



[(S z (t))] = S d (t) 



So(t) 



(3) 



where we have used [{S x (0))] r}lp = [r 2 ] r J2jti\Sz\j) = 
[r 2 ] r M T in the normalization; S d {t) = Mf 1 E,- Sjj(t) 



r Z^ji u 3i -J2 (t) are the diagonal and 



andS (t) = Mr'E^S; 

off-diagonal maxtrix element of electron spin. We note 
that Eq. ([3]) indicates that the averaged spin dynamics 
depends on the initial condition only via a single param- 
eter, (3 = ([r]l/[r 2 ] r ) \ [e~ lip ] v \ , which depends on the 
phase distribution function, V v , much more significantly 
than on the amplitude distribution function, V r , since 
[r] 2 ~ [r 2 ],. for the usual function of V r and r > 0. We 
note that although the experimental preparation of a co- 
herent nuclear spin bath (i.e. finite value of j3) is not 
easy at the present stage, it has been realized how to 
control the coherent electron spin dynamics via interac- 
tion with a single nuclear spin in diamond [111 ]. There- 
fore, at least in a small quantum dot system, a coherent 
preparation and control of a few nuclear spin can be still 
realized. In Fig. [3Jd) we show the time-evolution of both 
Sd(t) and S (t) of spin matrix element. Not surprisingly, 
they are of very different properties: the diagonal part 



[(S g (t))]* )V = M 2 [r 2 ] 2 S d {t) 2 according to Eq. © 
some algebra we can derive: [{S z (t)) 2 ] r v — [(S z (t))} 2 — 
([r% - 2[r 2 ] 2 r ) £ . 5 w (t) a + [r 2 ] 2 r £ 3 - , S^S U . Since we 
are interested in the upper bound of the NSD, we 
may use E^(*)<^(*) = E,-01& (t)P r S z (t)\j) < 

EiO'l&(t) 3 |j> ^ S ( 5 + 1 )E j 0'U) = S(S + 1)M T , where 
Pr = E; 10 CI i s to project a state onto the subspace T 
with electron spin eigenvalue S z — S. Here we have used 
the fact that for any state, the expectation value of S z (t) 2 
must be equal or smaller the expectation value of total 
electron spin S 2 , which is, however, a conserved quan- 
tity of our system (see Eq. (TT])). Similarly we also have 
E 3 S^(t) 2 < Eij \S id (t)f < S(5+l)M r , and therefore 
[0x(t)) 2 ]r, v -[0S))]r, v ~0(M r ). In other words, after 
devided by [(5,(0))]^ oc M 2 , we find a{t) cx Mp 1/2 and 
therefore goes to zero in the limit of N ^> 1 or Mr — > 00. 

We can also apply similar method to study the NSD 
of a phase coherent initial state, i.e. [e lip ] v / 0. 
After some algebra, the expansion of [(S z (t)) 2 ] ryV — 
[(S z (t))] 2 , will have two additional summations (be- 
sides of the two shown above) with nonuniver- 
sal prefactors: first, we have Ei,j,i = 
(V\S z (t)P r S z (t)\V) < (V\S z (t)S z (t)\V) < S(S + 
1)(V\V) - S(S + 1) = S{S + l)M r , where we 

define \V) = Ej \j) as an auxiliary state. Secondly, 

1/2 



M, 



(V\S z (t)P r S z (t)\V) < Ml ,2 J{V\S z (t) 2 \V) 



< 



\J S(S + I) Mr, where we have used the fact that the 
inner product of two vectors must be equal or smaller 
than the product of their length. Therefore, after renor- 
malized by the initial spin average, we find a(t) oc Mf 1 
and becomes to zero in large system size just as for the 
complete random phase case (£ = 0). From the above 
results, we conclude that no matter how much phase co- 
herence between spin eigenstates of the initial wavefunc- 
tion, the spin relaxation dynamics can be always univer- 
sal (with zero NSD) in the limit of infinitely large phase 
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FIG. 4: (a) Comparison between the calculated (S z (t)) (open 
circle) with £ = and 7 = 1 and {S z {t))oos (filled circle) 
obtained from density of states only, (b) Comparison between 
the calculated (S z (t)) with £ = 1 and 7 = 1 and the result 
obtained by including |-B max / m i n ) only (see text). 



space (Mr 3> 1). Similar derivation for dynamics of pure 
dephasing ({S x (t))) and Harm echo decay (p+_(r)) can 
be obtained straightforwardly. 



IV. MICROSCOPIC ORIGIN OF 
UNIVERSALITY 

After concluding the universality of the most general 
spin relaxation system (Eq. ([1])), in the rest of this paper 
we return to a less general case, zero magnetic field and 
nonintcracting spin bath (Hz = H n -n = 0), to study the 
microscopic origin of the universal spin relaxation curves 
shown in Fig. [^a)-(c). We first rewrite (S z (t)) in terms 
of energy eigenstates, \E): 



(Ss(t)) 



dEV(E) J dE'V{E')C E C* E , 
x{E\S z \E')e- l{E - E ")\ 



(4) 



where T>(E) is the density of states (DOS) of the system 
and Ce = (tpo\E). According to our numerical cal- 
culation, we observe that the matrix element, (E\S Z \E'), 
varies almost randomly for different energies. Since now 
we are interested in the simplest possible explana- 
tion for the features of electron spin dynamics, we 
can first neglect such structureless random ma- 
trix element for simplcity. As we will see later, 
it turns out that this simplication does bring a 
very useful understanding of the universal spin 
dynamics. 

For the case when initial wavefunction \ipo), is totally 
randomly distributed in the phase space, T, i.e. 7 = 1 
and £ = 0, we can further assume that (ipo\E) is also 
independent of energy E in above equation. As a re- 
sult, Eq. (fj| can be approximated by (S' z (i))DOS = 
\J dEV(E)e~ lEt \ , which is just power spectrum the den- 
sity of state. In Fig. 0Ja) we show the full numerical 
result of {S z (t)) for £ = compared to [S z (t))DOS given 
above. One can see that the later can qualitatively re- 
produce all the important structure of the full numerical 
results. This agreement helps us to conclude that the 
decay time of (S z (t)) is mainly determined by the width 
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FIG. 5: (Color online) (a) and (b) are the density of states 
of systems with different electron/nuclear spins. Note that 
in both figures, the energy axis for 1 = 1 cases (filled circle) 
have been rescaled by a factor 1/2 in order to fit the same 
scale as / = 1/2 case, (c) and (d) are ensemble-averaged 
spin relaxation curves for 5=1 and 7 = 1 case, with £ = 
and f = 1 respectively. For comparison, results of different 
nuclear spins are shown together after rescaling the time axis 
(see above). 



of DOS peaks (see Fig. [HJa)), while the time scale of the 
second peak of (S z (t)) is given by the energy separation 
between the two peaks in DOS. 

As for the single mode oscillation shown in Fig. [^c) 
for full spin coherent initial state (£ = 1), we can ap- 
ply similar study but notice that the coefficient, Ce, 
is not a constant for all eigenstate energy any longer. 
Our numerical calculation shows that \Ce\ is very small 

— 1/2 

(oc M r ) for all energy except near E = -E max / min , 



where \E, 



max/mm 



is the eigenstates of the top(maximum) 



and bottom(minimum) of the eigenenergy band. This is 
j) has a large overlap with some partic- 



because \E. 



max/mm 



ular spin eigenstate (as shown in Fig. [TJa) and (b)) and 
the overlapping coefficients are not cancelled out due to 
the same phase in a full spin coherent initial state = 1). 
In Fig. HJb), we compare the numerical result of the uni- 
versal dynamics (7 = 1) of a full coherent initial state 
(£ = 1) and the result calculated by using l-Emax/mm) 
only (with proper nomalization) . We find the agreement 
is excellent, predicting the same oscillation frequency and 
even the same phase. The agreement justifies the 
approximations used in the derivations after Eq. 
(4) and also shows that the universal behavior of the 
spin relaxation dynamics can be simply explained by the 
structure of density of states and the two special spin 
configurations as shown in Fig. QJa) and (b). As for re- 
sults with < £ < 1, it can be also explained well by a 
linear combination of above two results, as suggested by 
Eq. ©. 
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V. RESULTS FOR DIFFERENT SPINS 

After systematically investigating the spin relaxation 
dynamics for an spin-half electron inside a spin-half 
nucear spin bath, here we further extend the study of uni- 
versality to systems of different electron/nuclear spins. 
In Fig. |Sfa) and (b), we show the density of states 
for (5,7) = (S,I) = (±,1), (5,7) = and 

(S,I) = (1, 1) in different curves. For the convenience of 
comparison, we rescale the energy scale in each plot and 
normalize the hieght of DOS by the total size of phase 
space, T. Surpisingly we find that the DOS structure is 
almost the same for different nuclear spins I as long as 
the electron spin 5 is the same. This reflects the fact 
that the total Hilbert space of the nuclear spin bath has 
been large enough due to the number of nucluei so that 
the spin degrees of freedom does play very little role in 
the structure of energy spectrum. Analysing the energy 
eigenstate configuration, we find the spin configuration 
near the degeneracy regime (position of the peaks) are 
related to if the the central nuclei spin configuration is 
polarized and parallel (or anti-parallel) to the electron 
spin. Similar observation also applys to the triple peak 
structure in Fig. [5fb) for 5 = 1: In Fig. [5jc) and (d), 
we show the spin relaxation curves for 5 = 1 with spin 
phase random (£ = 0) and spin phase coherent (£ = 1) 
initial wavefunctions respectively, after properly rescal- 
ing the horizontal axis. One can see that results in (c) 
are very similar to the spin half case (Fig. Uta)), while 
it shows a beating oscillation for a coherent initial wave- 
function (d). The rescaled time-evolution for I = i and 
1 = 1 are very similar, except for a small phase twist. 
We then conclude that the the spin relaxation dynam- 
ics is insensitive to the nuclear spin degrees of freedom, 
consistent with our earlier statement that the universal 
spin dynamics is independent of the nuclear spin config- 
uration. Our results for 5 = I = 1 can be also applied to 
the study of spin dynamics in the mixtures of spinful cold 
atoms in all-optical trap, where the localized "electron" 
and the "nucluei" can be prepared easily by using optical 
lattice with proper wavelength difference. The advantage 
for cold atom system is that the initial spin configura- 
tion can be prepared easily and the coupling strength, 
Aq, can be tuned via optical Feshbach resonance and/or 
other method. 



VI. CONCLUSIONS 

In this paper we rigorously prove that the electron spin 
decoherence due to nuclear spin bath can be always uni- 
versal if only coupled by sufficient large amount of spin 
eignestates. There are several features about the uni- 
versal dynamics that we want to emphasize: First, in 
the derivation above, we do not rely on any particular 
form of the distribution function (V r and V v ), hence 
the universality of spin dynamics is independent of the 
nuclear spin configuration. However, if the initial state 



is composed of only finite numbers of spin eigenstates 
(as done in the literature), our derivation will fail since 
[r 2 ] r — > in the denominator of cr(t), i.e., no univer- 
sal dynamics can be expected. Secondly, the univer- 
sality does not rely on any particular Hamiltonian, so 
our conclusion also applies to systems which include nu- 
clear spin interaction, finite magnetic field, or any other 
more complicated system. Different system Hamiltoni- 
ans just bring different averaged results of spin dynam- 
ics, but the huge phase space (not necessary the huge nu- 
clear number) can always ensure it to be the most prob- 
able one regardless of the details in the initial condition. 
Such important results lead to another conclusion that 
a numerical simulation of a much smaller system (say 
N ~ 10 — 20) can still have large enough phase space 
(M n = (25 + 1) x (21 + 1) N ~ 10 3 ~ 6 ) and hence gives 
similar results as given by macroscopic number of nu- 
clei [3]. Excellent agreement between our results 
of small size calculation (Fig. [2](a) with N = 11) 
and a meanfield type calculation of a much larger 
system (for example, Fig. 4 of Ref. [8] with 
N = 2000) ensures the existence of such scale- 
indepedent universal dynsmics. Our results there- 
fore make a realistic comparison between a theoretical 
calculation and experimental data possible, leaving only 
a single unknown parameter, /3, as a fitting parameter. 
(For example, in the Fig. 4 of Ref. /3 = is 
expected due to the thermalized initial bath.) Fi- 
nally, our derivation relies on the fact that electron spin is 
a conserved quantity with an upper-bounded eigenvalue 
(not scaled with system size) . This may explain why spin 
eigenstate can be a special basis for studying universal 
physics and restricts a naive application of our results to 
the relaxation dynamics of other physical quantities. 

It is also worthy to note that the universal dy- 
namics may bear a close relationship with the 
quantum central limit theorem (QCLT). QCLT 
has been used to study the quantum state esti- 
mation without using a large ensemble [16] and to 
explain why quantum and classical random walks 
possess different behaviors [17]. It is nature to 
conjecture that the existence of the universal dy- 
namics and the reason why a small size system 
can already capture the behavior of the macro- 
scopic system can be understood in the context 
of QCLT. For example, in Ref [16], it was pointed 
out that a quantum state estimation with small 
error using small size ensemble is possible. This is 
clearly resemble to our work, where a small size 
system can capture the universal dynamics of a 
macroscopic system. However, such a connection 
is not at all transparent in the details of the ap- 
proach. Further study about such an interesting 
connection is beyound the scope of this paper but 
could be very interesting for furure investigation. 

We appreciate fruitful discussion with S. Das sarma, 
X. Hu, S. K. Saikin, L. J. Sham, and J. M. Taylor. This 
work is supported by NSC Taiwan. 
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